Dark Matter in B - L Extended MSSM Models 
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We analyze the dark matter problem in the context of supersymmetric U(1)b-l model. In this 
model, the lightest neutalino can be the B — L gaugino Zb-l or the extra Higgsinos xi,2 dominated. 
We compute the thermal relic abundance of these particles and show that, unlike the LSP in MSSM, 
they can account for the observed relic abundance with no conflict with other phenomenological 
constraints. The prospects for their direct detection, if they are part of our galactic halo, are also 
discussed. 



INTRODUCTION 

Nonvanishing neutrino masses and the existence of 
nonbaryonic dark matter (DM) are the most important 
evidences of new physics beyond the Standard Model 
(SM). A simple extension of the SM, based on the gauge 
group G B -L = SU(3)cxSU{2) L xU(l) Y xU(l)B-L,can 
account for current experimental results of light neutrino 
masses and their large mixing [1]. In addition, the extra- 
gauge boson and extra-Higgs predicted in this model have 
a rich phenomenology and can be detected at the LHC 
[2]. It is worth mentioning that several attempts have 
been proposed to extend the gauge symmetry of the SM 
via one or more U(l) gauge symmetries beyond the hy- 
percharge gauge symmetry [3, 4]. 

Within supersymmetric context, it was emphasized 
that the three relevant physics scales related to the super- 
symmetry, electroweak, and baryon minus lcpton (B — L) 
breakings are linked together and occur at the TeV scale 
[5]. Indeed, it was shown that radiative corrections may 
drive the squared mass of extra B — L Higgs from pos- 
itive initial values at the GUT scale to negative values 
at the TeV scale. In such a framework, the size of the 
B — L Higgs vacuum expectation value (VEV), responsi- 
ble for the B — L breaking, is determined by the size of 
the right-haneded Yukawa coupling and of the soft SUSY 
breaking terms. 

In this paper, we consider the scenario where the extra 
B — L neutralinos (three extra neutral fermions: U(1)b-l 
gaugino Zb-l and two extra Higgsinos x\,i) can be cold 
DM candidates. It turns out that the experimental mea- 
surements for the anomalous magnetic moment impose 
a lower bound of order 30 GeV on the mass of U(1) B -l 
gaugino Zb-l, while Higgsinos xi,2 can be very light. 
We examine the thermal relic abundance of these parti- 
cles and discuss the prospects for their direct detection 
if they form part of our galactic halo. 

It is worth mentioning that assuming the lightest 
neutralino in minimal supersymetric standard model 
(MSSM) as DM candidate implies sever constraints on 
the parameter space of this model. Indeed, in the case of 
universal soft-breaking terms, the MSSM is almost ruled 



out by combining the collider, astrophysics and rare de- 
cay constraints [6] . Therefore, it is important to explore 
very well motivated extensions of the MSSM, such as 
SUSY B — L model which provides new DM candidates 
that may account for the relic density with no conflict 
with other phenomenological constraints. 

The paper is organized as follows. In section 2 we 
briefly review the supersymmetric U(\)b-l model with 
a particular emphasis on its extended neutralino sector. 
Section 3 is devoted for computing the LSP annihilation 
cross section for Zb-l, Xi, and Xi- I n section 4 we exam- 
ine the possible constraints imposed by the experimental 
limits of muon anomalous magnetic moment on the mass 
of Zb-l- We discuss the relic abundance of these DM 
candidates in section 5. We show that they can account 
for the measured relic density without any conflict with 
other phenomenological constraints. The direct detec- 
tion rate of Zb-l and xi,2 is briefly discussed in section 
6. Finally we give our conclusions in section 7. 



U(1) B -l SUSY MODEL 

In B — L extension of the MSSM, the particle content 
includes the following fields in addition to the MSSM 
fields: three chiral right-handed superfields (iVj), a vec- 
tor superfield associated to U(1)b-l (Zb-l), and two 
chiral SM singlet Higgs superfields (xi, Xi)- This class 
of B — L extension of the SM can be obtained from a uni- 
fied gauge theory, like 50(10), with the following branch- 
ing rules for symmetry breaking: 50(10) is broken down 
to Pati-Salam gauge group: 5?7(4) c x SU{2) L x SU(2) R 
through the vacuum expectation value (vev) of the Higgs: 
(1, 1, 1) in 54# or 210jj representation at GUT scale. 
Then Pati-Salam can be directly broken down to B — L 
model: SU(3) C x SU{2) L x U(l) Y x U(1) B -l through 
the vev of the adjoint Higgs: (15, 1, 3) below GUT 
scale. Finally, the B — L model is broken down to 
5C/(3) C x SU{2) L x U(l)y at TeV scale as mentioned 
above. In this case, although the U(1)y and U(1)b-l 
are exact symmetries at high scale (larger than TeV), 
they are non-orthogonal. This can be seen by notic- 
ing that the orthogonality condition [7] is not satisfied: 
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J2f Y^ Yb-l 7^ 0, where Y? and Yg_ L arc the hyper- 
charge and the B — L charge of the fermion particle (/). 
In this respect, there is a kinetic mixing between the 
gauge fields of U(l)y and U(1)b-l- However, LEP re- 



sults [10] stringently constrain the corresponding mixing 
angle to be < 10~ 2 . Therefore, in our analysis, we ne- 
glect this small mixing and consider the following super- 
potential: 



W = {Yu)ijQiH u Uj + (Y D )ijQiHdDj + (Y L ) tJ L t H d E^ + (Y v )ijLiH u Nj + {Y N ) l3 N^N^ Xl + (J-(H u H d ) + p! X xXi- (1) 

I 



The B — L charges of superfields appeared in the super- 
potential W are given in Table I. 
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TABLE I: The [/(1)b-l charges of the superfields. 



For universal SUSY soft breaking terms at grand uni- by 
fication scale, Mx, the soft breaking Lagrangian is given 



J 



■-soft 



+ ^0 



iq.i 2 + m 2 + iai 2 + i^i 2 + m 2 + m 2 + \h u \ 2 + \H d \ 2 + \ X i\ 2 + \ X2 \ 

YuQU c H u + Y D QD c H d + Y E LE c H d + Y V LN C H U + Y N N c N c Xl 



1 



[B(fj,H u H d + n'xiXi) + h.c] + 2 M V2 



g a g a + W a W a + BB + Z b _ l Z B -l + h.c. 



(2) 



where the tilde denotes the scalar components of the chi- 
ral matter superfields and fermionic components of the 
vector superfields. The scalar components of the Higgs 
superfields H Utd and X i,2 are denoted as H Uid and xi,2, 
respectively. 

As shown in Ref. [5], both B—L and electroweak (EW) 
symmetries can be broken radiativcly in the supersym- 
metric theories. In this class of models, the EW, B — L 
and soft SUSY breaking are related and occur at the TeV 
scale. The conditions for the EW symmetry breaking are 
given by 



where 



ml + [i 2 , i = 1, 2 



m. 



-Bfi, 



< H u >= v u /V2, < H d >= v d /V2. 



tanp = — , 

Vd 

(4) 



Here rriH u and mn d are the SM-like Higgs masses at the 
EW scale. Mz is a neutral gauge boson in the SM. It is 
worth noting that the breaking SU (2)lxU (1)y occurs at 
the correct scale of the charged gauge boson (Mw ~ 80 
GeV). Similarly, the conditions for the B — L radiative 
symmetry breaking are given by [8] 



m 2 TT tan 2 8 



tan 2 8 - 1 



Mf/2, sin 2/? 



-2m% 



(3) 



z 2 [i\ — [i\ tan 2 9 
tan 2 6 - 1 



M| B _J2, sin 20 



-2/i§ 



(5) 
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where 

i = 1, 2 //§ = -B//, 
< Xi >= «i/>/2, < X2 >= v' 2 /V2. 



v' 

tan# = -4, 

«2 



gent constraint on U(1)b-l obtained from LEP 11 result, 
which implies [10] 



(6) 



M 



Zb- 



9B-L 



> 6 TeV. 



(7) 



Here m Xl and m X2 are the C/(l)s_L-like Higgs masses at 
the TeV scale. The key point for implementing the radia- 
tive B — L symmetry breaking is that the scalar potential 
for xi an d Xi receives substantial radiative corrections. 
In particular, a negative squared mass would trigger the 
B — L symmetry breaking of£/(l)s_£. It was shown that 
the masses of Higgs singlets xi an d Xi run differently in 
the way that m Xl can be negative whereas m 2 2 remains 
positive. The renormalization group equation (RGE) for 
the B — L couplings and mass parameters can be de- 
rived from the general results for SUSY RGEs of Rcf. 
[9]. After B — L symmetry breaking, the U(1)b-l gauge 
boson acquires a mass [1]: M\ = Ag 2 B _ L v' . The high 
energy experimental searches for an extra neutral gauge 
boson impose lower bounds on this mass. The most strin- 



Now we analyze mass-spectrums which have some de- 
viations from MSSM-spectrums in particular, SM singlet 
Higgs bosons, the right-handed sncutrinos, and the neu- 
tralinos. The Higgs sector in the SUSY B — L exten- 
sion of the SM consists of two Higgs doublets and two 
Higgs singlets with no mixing. However, after the B — L 
symmetry breaking, one of the four degrees of freedom 
contained in the two complex singlet xi an d Xi are swal- 
lowed by the Z B _ L to become massive. Therefore, in 
addition to the usual five MSSM Higgs bosons: neutral 
pscudoscalar Higgs bosons A, two neutral scalars h and 
H and a charged Higgs boson H^, three new physical 
degrees of freedom remain [5] . They form a neutral pseu- 
doscalar Higgs boson A' and two neutral scalars h! and 
H' . Their masses at tree level are given by 



J 



rn 2 H >, h , = \ {ml + M% B _ L ± 



) 2 -4m^M| B [ cos26») . 



(8) 



The physical CP-even extra-Higgs bosons are obtained 
from the rotation of angle a: 



h' 
W 



cos a sin a 
— sin a cos a 



Xi 

X2 



where the mixing angle a is given by 



1 

a = - tan 
2 



M'? 

tan 29 



M'p 



M'l 



M' 2 



(9) 



(10) 



For v[ » v' 2 , one finds the mixing angle a is very small, 
hence the above diagonalizing matrix is close to the iden- 
tity. In this case, to a good approximation, one can as- 
sume that hi = xi an d H' = X2- We are going to adopt 
this assumption here. 



Now we turn to the right-handed sncutrinos, in the 



basis of {(j>„ 



4>v 



with (j> v 



(v L , v* L ) and 



{vr, i>x), the sneutrino mass matrix is given by the fol- 
lowing 12 x 12 hermitian matrix: 



2 _ I Ml LUL Ml LVR 
J ~ M 2 M 2 

" ± VBML VRVR 



where M 2 AVB (A, B = L,R) can be written as 



(11) 



Mi 



M %B M A*B 
M %B M AfB 



with 



(12) 



J 



M\ 
M 2 t 

V R V R 

M 2 T 
M 2 T 

V R VR 

M 2 t 



M 2 



2/3 + v 2 sin 2 l3Ui INS (YjY„)U M NS 



U MNS m U MNS + ~Y~ C0S 

m 2 + M 2 Nl 

v sin (3U f MNS A (Y N y + v cos, ft idJ ] MNS A Q (Y v ) ] , 

v'sm.6A {Y N )\ 

wsin A (Y V )M N , 



M t = 0, 



4 

(13) 



where v' sm.9 =< \i >, Mjy — Y^v' and Umns is 3 x 3 
unitary matrix termed the Maki-Nakagawa-Sakatalepton 
mixing matrix [12]. Therefore, in general the order of 
magnitude of the sneutrino mass matrix is as follows: 



M 



0(v 2 ) 0(vv') 
0(vv') 0{v' 2 ) 



(14) 



Since 



TeV, the sneutrino matrix elements arc of the 



same order and there is no seesaw type behavior as usu- 
ally found in MSSM extended with heavy right-handed 
neutrinos. Therefore a significant mixing among the 
left- and right- handed sneutrinos is obtained. The phe- 



nomenological consequences for such mixing have been 
studied in [13]. 

Finally, we consider the neutralino sector. The neutral 
gaugino-higgsino mass matrix can be written as: 



M 7 (B, W 3 , H° d , Hi xi, X*, Z B -l) 



M A O 
O t M 3 J' 

(15) 



where the AI4 is the MSSM-type neutralino mass matrix 
and M 3 is 3 x 3 additional neutralino mass matrix, which 
is given by: 





M,= \ -fx' 

—Igs-hv' sin# 



—/j,' — 2gB-Lv' sin 6 

2gs-Lv' cos9 

2g B -Lv'cos9 M 1/2 



(16) 



As a feature of the orthogonality of U(l)y and U(1)b-l 
in this class of models, there is no mixing between M4 
and At 3 at tree level. Note that in extra U(l) gauged 
models, which are proposed to provide an explanation 
for the TeV scale of /Lt-term through the vev of a sin- 
glet scalar, the neutralino mass matrix is given by 6 x 6 
matrix. If the the extra singlet fermion is the lightest 
neutralino, then it can be an interesting candidate for 
dark matter, as shown in Ref.[14]. In our case, one di- 
agonalizes the real matrix M7 with a symmetric mixing 
matrix V such as 



other coefficients are close to zero. In our analysis, we 
will focus on these two types of LSP and analyze their 
potential contributions to DM in the universe. The 
mass eigenstates of the matrix M3 are in general non- 
trivial mixtures of the fermions (xi, X2, Zb-l)- The 
limit of pure Zb-l that we consider can be obtained if 
v' « fj,' and the limit of pure Xi(2) can be obtained if 
fi', v' sin(cos)# << v' cos(sin)^ 1 . 



VMjV T — diag.(m 0), k = l,..,7. 



(17) 



In this aspect, the lightest neutralino (LSP) has the fol- 
lowing decomposition 

X? - V n B+V 12 W 3 +V 13 H° d +V u H°+V 15 xi+V 16 X2+V 17 Z B ^ 

(18) 

The LSP is called pure Zb-l if Vn ~ 1 and Vu <~ 0, 
i = 1,..,6 and pure Xi(2) if ^15(6) ~ 1 an( i au the 

1 We would like to thank the referee for drawing our attention to 
this point. 



LSP ANNIHILATION CROSS SECTION IN 
U(1) B -l SUSY MODEL 

L • 

As advocated in the previous section, we focus on the 
cases where LSP is pure Zb-l or Xi(2)- La this case, the 
relevant Lagrangian is given by 



-jC~ b l ~ tV2g B - L Yl_ L Z B -LPR.ffL + iV2g B -LYl_ L Z B -LPLffR + ex., (19) 
-£~ ~ z^. 9B _ L r^ L ^ 1 ^_ L7 5Xi+^.9B-L^i L Xi^-LXi + (>V) 4J Xi^ c ^; + cc, (20) 
-C~ ~ iy/2g B -LY]?_3 2 pB-Lri)b + iV2g B -LY^_ L %Z B -LX2 + cc, (21) 
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where / refers to all the SM fermions, including the 
right-handed neutrinos, /l and are the left-handed 
and right-handed sfermions mass eignstates respectively. 
Yg_ L is the B — L charge defined in Table I. We assume 
the first right-handed neutrino N% is of order 0(100) 
GeV, therefore the annihilation channel of the LSP into 
NiNi is also considered. 



■*■ x., 



- ► x 




FIG. l^The dominant annihilation cross sections in the case 
of the .ZB_L-like LSP. Note that u-channel is also taken into 
consideration for each of diagram. 



From Eq. (19), one finds that the dominant annihila- 
tion processes of xl = Zb-l are given in Figure 1. Our 
computation for the annihilation cross section leads to 
the following a 



tion (<7 n 



z B - 
a 



and b~ , where the approxima- 
tes-!, 

b v 2 , with v is the velocity of the 



incoming LSP, is assumed: 



Zb-i 



Zb-i 



9b-l 
547rm 2 



[(3' t r 2 z 2 + 27/3' N r 2 N z%] 



(22) 



1627rm o 



(1 - 2r f + 2rf) + 



9b-l 
Aitm 2 a 



r P'r 2 

1 {ai + r\ + z 2 (a 4 + r 4 )} t + t r\{a\ + r\ + z 2 (a 4 + r 4 )} 



27 



N 



^^\O ln Of m \ 2 \V 2 , l+4 V 2 T l+4 \ 2 (3' x r 2 x 

x 2 -3 



Xi 



;(l + w 2 )p: 



■12 



(i = 1 or 2) 



(23) 



ai = 



2 2^2 _ __ 

3 + 4 x a z a i2 z a , a*- 4 



ri = Y [-4 + z 2 a + r a (4 - 'iz 2 a - z A a )] , r 4 = ^(-3 + 2z 2 a + 5r a (3' 2 ), 



m a /m x o, w a = m a /m x o, r a = (1 - z a + w a ) , a = 1 - z a , x a 



2(1-^)' 



(24) 



where m a is a final-state mass, m a is a mediated-particle 
mass, O is the extra Higgs mixing matrix, as defined 
in Eq. (9) and Eq. (10). In our approximation, Oi n 
is given by Oi n — §i n , and we set m j = m^ R ~ m f L ■ 
Moreover, V 2i is the coefficient of Next LSP (NLSP). We 
assume that Xi(2) is our NLSP, therefore V 2 ,i+4 — l(i 
I or 2), V 2 j ~ 0(j : ^ 5 or 6). In the range of parameter 
space that we consider, the values of a— and b~ 

Zb-l Z B -l 

are typically < 10~ 8 . For m~ > 100 GeV, the an- 



nihilation channels into extra-Higgs xi an d Xi may give 
the dominant contributions to b~ 

Zb-l 

Now we turn to the Higgsino contributions. From Eq. 

(20), one finds that the dominant annihilation processes 

of Xi = Xi are given in Figure 2. The computation of 

the cross section leads to the following results for a~ and 

to Xl 

b~ : 

Xl 



P'n z n 



9%-L m x o 



^/2{Yn,1t- 



)Wb- 



4m| B 



m 2 

Xi 



l Z B - 



(25) 



6~ = 



3tt 



7T 



0W)V 



N 



x? 



(ai + n)iv + 



xi 



1 - 4 



aiAr 



FIG. 2: The dominant annihilation channels of the Xi-like 
LSP. For the last two diagrams, the it-channel is also consid- 
ered. 



ft 



N 



V2(Y NAm ) 2 g%_ L 



I '7 



( i 2 z qI2\ 

{ax + r 5 z - -rp )n 



PmOl 



\ 2 \V 2 rV^P' x rl 



1" 



(26) 



Here Vn is the coefficient of NLSP. We assume that Zb-l 
is our NLSP, therefore V 27 V 2i ~ 0(« ^ 7). We also 
assume (ljv)i m ~ (3^v)n. 



FIG. 3: The dominant annihilation cross section in case of 
the X2-like LSP. Note that u-channel is also taken into con- 
sideration for the t-channel diagram. 



Zb-l 



Finally we consider the annihilation process of \i = 
Xi- From Eq. (21), one finds that X2X2 annihilation is 
dominated by the diagrams in Figure 3. The computa- 
tion to the cross section of \2 leads to = and b 
is given by 



V2 



X'2 



b~ 

X2 



3tt 

4 



Mi 

Z<B~L 



3 2 V tJ 



_ T 



\O2nOi 



2 

x' x 



1 + - x 2 z 2 - —z 2 
3 4 * * 12 * 



rrr 



1-4 



(27) 



It is remarkable that for m~ > 100 GeV, their anni- 
hilations are dominated by extra-Higgs channel. There- 
fore, b~ is very close to b~ and a~ is quite suppressed. 

X 1 . . X2 XI _ 

Thus, in this region of parameter space both \i an( i Xi 
have very similar annihilation cross section values. 



CONSTRAINTS FROM MUON ANOMALOUS 
MAGNETIC MOMENT 

In the case of Ze_£-like LSP, a significant contribu- 
tion to muon anomalous magnetic moment (a^) may be 
obtained due to the 1-loop diagram mediated by Zb-l 



and smuon, as shown in Figure 4. Note that X1.2 Lave 
no direct couplings with the SM fermions, thus they do 
not contribute to a M . The recent experimental value has 
been determined with a very high precision by the E821 
Collaboration at the National laoratory [16] 

a^ p - = (116592080 ± 60) x KT U . (28) 

This value differs from the SM predicion by the following: 

Aa p = a e * p - - al M = (278 ± 82) x lO" 11 . (29) 

Therefore, Zb-l contribution to a M should satisfy the 
following constraints: 

1.96 x 10~ 9 < Aa^ B - L < 3.6 x 10" 9 . (30) 




Zb-l 



Our computation for Zb-l contribution to a M leads to 
the following result: 



FIG. 4: Zb-l contribution to the muon anomalous magnetic 
dipole moment 



Aaf, B - L 



9 2 B-L 



E 



7T 2 ^ 6mf (1 - Si y 

i=l,2 Mi v % > 



x - Si )(C/ A ) 2i (C/ A ) li Yi_ i Y#_ i 6(l - ,s 2 + 2* In*) 



+ -^(lOWlUl 2 + I(^)i^A-l| 2 )(1 - 6 Sl + 3s 2 + 2s\ - 6*? In*)], 



(31) 



where [7^ is a diagonalized unitary matrix of the slep- 
ton sector, s* = (m x a /m mji ) 2 and Y B _ L is the ?7(1)b_l 
charge in the Table I. This result is consistent with the 
derivation of the new contribution to in supersymmet- 
ric U(iy model [17]. 

Here few comments are in order: (i) The second term 

in Aafi, B ~ L is suppressed by m l _ t /mp i ~ 0(1O -3 ), while 
the first term is proportional to the off-diagonal elements 
of the diagonalized matrix Up which are typically of order 
0(1(T 2 ). Therefore the first term is Eq. (31) gives the 

dominant contribution to Aa M . (ii) From the Eq. (30), FIG. 5: Aa^ B - L versus the mass of the LSP Z B -l for mp, R 
the sign of Zb-l contribution to Aa^ should be positive. 
Thus [(Up)u(Up) 2 i + (Up) 12 (Up) 2 2]Y B _ L Y£_ L must be 
negative. Note that Si < 1, hence the function f(si) = 
\/Sj(l — s i)(l — s f + 2si In Si) is always positive. The ele- 
ments of Up have a sign difference that helps in satisfying 
this requirement and allows for positive contribution to 




20 40 60 80 100 120 140 160 180 200 



m.(GeV) 



m jlL ~ A ~ 300 GeV, /x ~ 500 GeV with tan/3 = 10, 20 
and 30, and Qb-l ~ 0.5. 

< m x o < 100 GeV. While for tan/3 = 20, the allowed 
region of m^o is rather wider: m^o > 60 GeV. 

° Xi Xi r*j 



Ac 



For example, in case mpL = mp.R = A 



300 



GeV, \i — 500 GeV and tan/3 = 10, the corresponding 
Up matrix is given by 



LSP RELIC ABUNDANCE IN U(1) B -l SUSY 
MODEL 



2 1 



2 1 



(32) 



(iii) Large values of tan/3 enhance the off-diagonal ele- 
ments of Up. Hence Aaf L B ~ L are enhanced by large values 
of tan/3. 

In Figure 5, we plot Aafi B ~ L as a function of the Zb-l- 



like LSP mass, 



>, for tan/3 = 10, 20 and 30. Other 



SUSY parameters are fixed as above. From this figure, it 
can be easily seen that a significant B — L contribution 
to Aa M can be obtained for tan/3 > 10. For tan/3 = 30, 
the LSP mass is constrained within the region 30 GeV 



In this section, we compute the LSP relic abundance in 
U(1)b-l SUSY Model. We adopt the standard compu- 
tation of the cosmological abundance, where the LSP is 
assumed to be in thermal equilibrium with the SM par- 
ticles in the early universe and decoupled when it was 
non-relativistic. Therefore, the LSP density can be ob- 
tained by solving the Boltzmann equation [15]: 



dn 



dt 

where 



Xl - + 3Hn x « = -< a a x Tv > [(n x? ) 2 



(33) 



is LSP number sensity with m^o = a v ori v o. 
One defines f2 x o = p x o / p c , where p c is the critical mass 
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FIG. 7: Q.h 2 versus Xi, 2 -Uke LSP mass for g B -L G [0.1,0.5] 
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FIG. 6: Qh 2 versus Z B -L-Hke LSP mass for g B ~L S [0.1, 0.5]. 

density. It turns out that [15] 

3.76 x IQ-^GeV- 2 



xp = In 



0.0955m p ;m x o i (a + 6b/xp) 

(g*(T F )x F ) 1/2 



(34) 



where m p i is the Planck mass (1.22 x 10 19 GeV) and 
9*(T F ) enumerates the degrees of freedom of relativistic 



particles at Tp. From the expressions, one notes that the 
LSP relic abundance depends only on the LSP mass and 
the annihilation cross section coefficients a and b. 

In our numerical calculation for the LSP annihilation 
cross section, we consider the following values of masses 
for the particles contributing in the process (extra-light 
Higgses (x?( 2 ))> s f erm i° ns (/)> the lightest right-handed 
neutrino (Ni) and the NLSP (x°)): m Xl(2) = 100 GeV, 
/ = 200 GeV, Ni = 100 GeV, m x » = m x? + 30 GeV. 

In Figure 6, we present the values of relic density 
Q.h 2 as a function of the LSP mass for Zs-L-like LSP 
and gs-L € [0.1,0.5]. The horizontal lines are experi- 
mentally allowed regions from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [18] results for cold dark 
matter relic density. Here, we have imposed the con- 
straint on the mass of Zs_L-hke LSP: m~ > 30 GeV, 
due to the experimental limits on the muon anomalous 
magnetic moment. From this figure, one notes that since 
the sfermion mass is fixed at 200 GeV, the annihila- 
tion channel due to its exchange produces a resonance 
at the LSP mass of order 100 GeV. For g B - L = 0.2, 
the allowed region is rather wide: 130 GeV < m x o, 
while for gs-L = 0.3, the allowed region is reduced to 
around 120 GeV. Finally, for g B - L > 0.4 a lighter LSP 
( m z ~ 100 GeV) is favored. 



[htb] 



Now we turn to the Higgsino X1.2 LSP. In Figure 7, we 
plot the LSP relic density £lh 2 as a function of \i or X2 
mass. As expected the relic abundance of xi or X2 are 
quite similar since they have very close annihilation cross 
section. From this figure, one notes that for gs-L < 0.1, 
there is no essentially any allowed region due to the fact 



that the relic abundance becomes quite large. While for 
gs-L = 0.2, 0.3, 0.4, 0.5, the allowed regions for m Xl 2 are 
given by [150 - 190] GeV, [130 - 135] GeV, [125 - 130] 
GeV, and [115- 120] GeV, respectively. 
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LSP DETECTION RATE IN U(1) B -l SUSY 
MODEL 



In this section we analyze the effect of the event rates 
of our relic neutralinos (Zb-l, Xi(2)) scattering off nuclei 
in terrestrial detectors. The direct detection experiments 
provide the most natural way of searching for the neu- 
tralino dark matters. The differential cross section rate 
is given by [19] 



dR 
dQ 



2m^om 2 . 

Xi r 



\F(Q)f 



dv, 



(35) 



where/i(i>) is the distribution of speeds relative to the 

m x am N 

detector. The reduced mass is m r = . where 



rriN is the mass of the nucleus, v n 



m +m N ' 
1/2 



[ Qm N \ 

\ 2m r J 



the energy deposited in the detector, and p x < 
sity of the neutralino near the Earth. It is common to 
fix p x o to be p x o — O.iGeV/cm 3 . The quantity a is the 



elastic-scattering cross section of the LSP with a given 
nucleus. In our model, a has two contributions: spin- 
independent (scalar) contribution due to the squark ex- 
change diagrams for Zs-L-^e LSP, and spin-dependent 
contribution arising from Zb-l gauge boson exchange di- 
agrams for x 1 ( 2 )-like LSP. For 76 Ge detector, where the 
total spin of 76 Ge is equal to zero, we have a contribution 
from the scalar part only, which is given by 



Ami 



\Zf p + (A-Z)f n 



(36) 



where Z is the nuclear charge, and A — Z is the number 
of neutrons. The expressions for the effective couplings 
to proton and neutron, f p and /„, can be found in Ref. 
[20]. Finally, the form factor F(Q), in this case is given 
by [19] 



F b (Q) 



3ji(gfli) 
qRi 



-q s 



(37) 



where q = ^/2m^Q is the momentum transferred and R\ 
is given by Ri = (R 2 - 5s 2 ) 1 / 2 with R = 1.2/mi 1 / 2 and 
A is the mass number of 76 Ge. ji is the spherical Bessel 
function and s ~ lfm. 

For 73 Ge detector, where the total spin of 73 Ge is equal 



to J 



we have a contribution from spin-dependent 



part only, which can be written as 



(Q)\ 2 = 



Am 



2 J 



'-[(fpYSM + U a n YS nn {q) + fpfnSpniq)], 



(38) 



where S pp (q) = S Q0 (q) + S n (q) + S 01 (q), S nn (q) = 
Soo(q) + Sn(q)-S i(q) and S pn (q) = 2[S o(q)- Six], and 
the expressions for /" and /" can be found in Ref. [21]. 



(events/kg/day) 

1E-5 



9=. L = - 5 
9b-l= - 4 



9b- l = - 3 



g =0.1 

B-L 



10 100 120 140 160 

Z-ino like LSP mass(GeV) 



is FIG. 8: Detection rate versus Zs-L-like LSP mass for gs-L 6 



is the den- [0.1, 0.5]. As in previous figures rriq — 200 GeV is assumed. 



The values of the spin structure functions Sao(q), Sn(q) 
and 5*oi (g) are given in [20]. 

In case of i?B_£-like LSP, the effective couplings to 
proton and neutron are very similar i.e. f p ~ /„. There- 



fore, the cross section, a SI = a— 1 , is given by 



t5 j 

Zb- 



Am l r 

7T 



^ 1 <iV|^|iV>^^ 



.xq9qn k 



^ 2 

i 



k=l 



1h 



. (39) 

where q refers to u, d, s, c, b, t. The hadronic matrix ele- 
ments are given by < N\qq\N >= f^m p /m q . The values 
of the parameters can be found in Ref. [21]. From 

Eq. (19), one finds that Z B -l couples universary to all 
type of quarks, i.e. 9q LhXq = 9q RkXq ^ iV2g B -LY^_ L . 

In Figure 8, we present our numerical results for the 
event rate R as a function of ^s_L-hke LSP mass for 
rn q = 200 GeV and gB-L G [0.1,0.5]. As can been seen 
from this figure, the detection rates are quite sensitive to 
the value of gauge coupling gB-L- This is due to the fact 
that R depends on the forth power of gB-L- Neverthe- 
less, the detection rates are less than 10 -3 events/kg/day, 
which are below the current experimental limit: 0.01 
events/kg/day [24] . Thus, one can conclude that Zb-l 
is beyond the reach of near future experiments. 

Now we turn to the case of Xi(2)-like LSP. As men- 
tioned above, in this case the scattering cross section is 
given by the spin-dependent part: <r SD = a~ D , which 

Xl(2) 

is given by Eq. (38) with 
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q=u,d,s \ ^B-L / \ / q=u,d,s 



Here we have used the lower limit on the ratio: 
Mzb-l I 9b-l reported in Eq. (7). The numerical values 
of 500(9), Sn(q), Soi(q) and (Aq)iv can be found in Ref. 
[21]. From this expression, it is clear that the detection 
rates of the extra Higgsinos-like LSP are extremely small. 
They are typically less than 10~ 16 (events/kg/day). This 
result is consistent with the spin-dependent contribution 
for the singlino in SUSY models with [7(1)' [22, 23]. How- 
ever, in this class of model, unlike our U{V)b-l model, 
the singlino dominated LSP may imply large detection 
rates, due to the spin independent contributions. 

CONCLUSIONS 

We have studied the DM problem in supersymmct- 
ric B — L extension of the SM. We showed that the 
extra B — L ncutralinos (three extra neutral fermions: 
U(1)b—l gaugino Zb-l and two Higgsinos Xi.,2) are in- 
teresting candidates for cold DM. We provided analytic 
expressions for their annihilation cross sections. Wc also 
computed the Zb-l contribution to muon anomolous 
magnetic moment and showed that the current experi- 
mental limits impose a lower bound of order 30 GeV on 
Zb-l mass. We analyzed the thermal relic abundance of 
both Zb-l and xi,2- We showed that unlike the LSP in 
MSSM, these particles can account for the measured relic 
abundance with no conflict with other phcnomenological 
constraints. Finally, we discussed their direct detection 
rates and showed that they are beyond the reach of our 
near future experiments. 

Acknowledgments 

We thank D. A. Dcmir for discussions. H.O. would like 
to thank Y. Daikoku, S. Nakamura and Y. Kajiyama for 
interesting comments. This work was partially supported 
by the ICTP grant Proj-30 and the Egyptian Academy 
for Scientific Research and Technology. 



[1] S. Khalil, J. Phys. G 35, 055001 (2008) [arXiv:hep- 
ph/0611205]. 

[2] M. Abbas and S. Khalil, JHEP 0804, 056 (2008) 
[arXiv:0707.0841 [hep-ph]]; W. Emam and S. Khalil, 
Eur. Phys. J. C 52, 625 (2007) [arXiv:0704.1395 [hep- 
ph]]; K. Huitu, S. Khalil, H. Okada and S. K. Rai, 



r 



arXiv:0803.2799 [hep-ph], to appear at Phys. Rev. Lett.; 
S. Khalil, O. Soto, arXiv:0804.0336 [hep-ph], to appear 
in JCAP. 

[3] Vernon D. Barger, Ernest Ma, and K. Whisnant, Phys. 
Rev. D 26, 2378(1982). 

[4] R. N. Mohapatra and R. E. Marshak, Phys. Rev. Lett. 
44, 1316 (1980); R. E. Marshak and R. N. Mohapa- 
tra, Phys. Lett. B 91, 222 (1980); C. Wetterich, Nucl. 
Phys. B 187, 343 (1981); A. Masiero, J. F. Nieves and 
T. Yanagida, Phys. Lett. B 116, 11 (1982); R. N. Moha- 
patra and G. Senjanovic, Phys. Rev. D 27, 254 (1983); 
R. E. Marshak and R. N. Mohapatra, In *Gotsman, E. ( 
Ed.), Tauber, G. ( Ed.): From Su( 3) To Gravity*, 173- 
181; W. Buchmuller, C. Greub and P. Minkowski, Phys. 
Lett. B 267, 395 (1991). 

[5] S. Khalil, A. Masiero, Phys. Lett. B 665: 374-377, 
(2008). 

[6] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, 
Phys. Lett. B 565, 176 (2003) [arXiv:hep-ph/0303043]. 

[7] Will Loinaz, Tatsu Takeuchi, Phys. Rev. D60 (1999) 
115008 [hep-ph/9903362 ]. 

[8] S. Khalil, Q. Shan, Nucl. Phys. B 564: 19-28, (1999). 

[9] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50, 2282 
(1994)[arXiv: hep-ph/9311340]. See also K. Inoue, A. 
Kakuto, H. Komatsu and S. Takeshita, Prog. Theor. 
Phys. 68, 927 (1982)[Erratum-ibid. 70, 330 (1983)]; J. 
R. Ellis, J. S. Hagelin, D. V. Nanopoulos and K. Tam- 
vakis, Phys. Lett. B 125, 275 (1983); L. E. Ibanez and C. 
Lopez, Phys. Lett. B 126, 54 (1983); L. Alvarez-Gaume, 
J. Polchinski and M. B. Wise, Nucl. Phys. B 221 495 
(1983); T.Kikuchi, T. Kubo, Phys. Lett. B 666: 262-268 
(2008). 

[10] M. Carena, A. Daleo, B. A. Dobrescu and T. M. P. Tait, 

Phys. Rev. D 70, 093009(2004) 
[11] F. del Aguila, M. Masip and M. Perez- Victoria, Nucl. 

Phys. B 456, 531 (1995); M. x. Luo and Y. Xiao, Phys. 

Lett. B 555, 279 (2003); P. H. Chankowski, S. Pokorski 

and J. Wagner, Eur. Phys. J. C 47, 187 (2006). 
[12] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 

28, 870 (1962) 

[13] Y. Grossman and H. E. Haber, Phys. Rev. Lett. 78, 3438 

(1997) [arXiv:hep-ph/9702421]. 
[14] D. Jarecka, J. Kalinowski, SF. King, and J. P. 

Roberts, ECONFC0705302: SUS15, (2007) arX- 

ive:0709.1862; J. Kalinowski, SF. King, and J. P. 

Roberts, arXiv:0811.2204. 
[15] K. Griest, Phys. Rev. D 38, 2357 (1988). 
[16] G. W. Bennett et al, [Muon G-2 Collaboration], Phys. 

Rev. D 73, 072003 (2006) [arXiv:hep-ph/9511378]; 

arXiv:hep-ph/9707451. 
[17] Y. Daikoku, hep-ph/0107305. 

[18] C. L. Bennett et al, [WMAP Collaboration], Astrophys. 
J. Sippl. 148, 1 (2003); D. N. Spergei et al, [WMAP 
Collaboration], Astrophys. J. Sippl. 148, 175 (2003) 

[19] S. Khalil, A. Masiero and Q. Shan, Phys.Rev. D 56: 
5754-5760,1997. 



11 



[20] Gerard Jungman, Marc Kamionkowski and Kim Griest, 

Phys.Rept. 267:195-373,1996. 
[21] P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, Mia Schelke 

and E.A. Baltz, JCAP 0407: 008, 2004. 
[22] V. Barger, P. Langacker, lan Lewis, Mat McCaskey, 

G. Shaughnessy and B. Yencho, Phys.Rev. D75: 



115002,2007. 

[23] B. de Carlos and J. R. Espinosa, Phys. Lett. B 407, 12 

(1997) [arXiv:hep-ph/9705315]. 
[24] Z. Ahmed et al. [CDMS Collaboration], arXiv:0802.3530 

[astro- ph]. 



